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Abstract 

We address the scaling behavior of the covariance of the magnetic field in the 
three-dimensional kinematic dynamo problem when the boundary conditions 
and/or the external forcing are not isotropic. The velocity field is gaussian 
and (^-correlated in time, and its structure function scales with a positive 
exponent The covariance of the magnetic field is naturally computed as 
a sum of contributions proportional to the irreducible representations of the 
S0(3) symmetry group. The amplitudes are non-universal, determined by 
boundary conditions. The scaling exponents are universal, forming a discrete, 
strictly increasing spectrum indexed by the sectors of the symmetry group. 
When the initial mean magnetic field is zero, no dynamo effect is found, 
irrespective of the anisotropy of the forcing. The rate of isotropization with 
decreasing scales is fully understood from these results. 
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I. INTRODUCTION 



The aims of this paper are two- fold. First, we are interested in the statistical properties 
of magnetic fields advected by turbulent velocity fields. Such magnetic fields possess a "self- 
stretching" term that is absent in the context of advected passive scalars (for a general 



and anomalous scaling needs to be addressed. Secondly, we want to focus on the anisotropic 
nature of turbulence: generically turbulence is forced by agents that are neither isotropic nor 
homogeneous, but most of the fundamental theories regarding universal scaling properties 
consider an ideal model of isotropic turbulence. In the case of a magnetic field advected 
by a Gaussian, 5-correlated velocity field with nontrivial spatial scaling we can present an 
exact (nonperturbative) solution of the full spectrum of anomalous scaling exponents of all 
the anisotropic contributions to the covariance of the magnetic field. We can thus offer a 
precise picture of the rate of isotropization upon diminishing scales, assess the importance 
of anisotropy for "inertial range" scaling, etc. 

The equation of motion of a magnetic field B{r,t) reads 



where u is the advecting velocity field, and / is the external forcing, and k is the mag- 
netic diffusivity. We address a model in which the velocity is taken Gaussian, isotropic, 
5-correlated in time, and its correlation function is 




Thus a dynamo effect may exist, and its relation to intermittency 



dtB{r, t) + u{r, t) ■ VB{r, t) = B{r, t) ■ Vn(r, t) 





(M"(r, t)u^{r', t')) = 6{t - t')D''^{r - r') 





The structure function S scales with exponent ^, < ,^ < 2: 




A < r < A . 
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On the other hand the forcing / is taken here Gaussian, (5-correlated in time but non- 
isotropic. The correlation function of the forcing has compact support in fc-space in an 
interval < k < 1/L, where L is the outer scale of the forcing /. We denote F°'^{R) = 

We are interested in the properties of the covariance of B, C°'^{R,t), 

C'^^iR,t)^{B'^iR,t)B^iO,t)) , (1.4) 

and eventually in the stationary quantity C^^{R) which is obtained in the stationary state if 
the forcing is balanced by dissipation. The calculation of this object in an isotropic ensemble 
has been presented by Vergassola 0. The anisotropic problem has been addressed recently 
by Lanotte and Mazzino 0. In the latter study the covariance ( |1.4|) was not properly 
expanded in terms of irreducible representation of the SO (3) symmetry group, and therefore 
an apparent mixing of the different sectors was found. As a result the authors had to tackle 
an infinite set of equations for all the sectors of the symmetry group. We show below that 
this mix-up is spurious, originating from an improper expansion. Due to the improper choice 
of the expanding basis set, an infinite linear coupling between all sectors appeared in 0. 
In order to solve the infinite linear system the authors were forced to assume the existence 
of a hierarchy between exponents belonging to different sectors and then only a posteriori 
to check the correctness of their assumption. In this way the calculation ends up with one 
correct set of exponents, as is shown below by using the proper expansion. We compute 
additional exponents that were left out in . We will also concern ourselves with the issues 
of the dynamo effect and the attainment of a stationary solution for (|1.4| ). 

The structure of this paper is as follows: in Sect. 2, after presenting the equations of 
motion of the covariance, we expand the solutions in terms of basis functions of the SO (3) 
symmetry group. In Sect. 3 the above expansion is used to obtain the matrix representation 
of the linear operator which determines the dynamics of the covariance. In Sect. 4 we use 
this matrix representation to show the absence of a dynamo effect in the anisotropic sectors 
of the covariance. Sect. 5 is devoted to the calculation of the anomalous scaling exponents 
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in the anisotropic sectors, and Sect. 6 offers a summary and a discussion. 



II. BASIC EQUATIONS AND THE DECOMPOSITION IN TERMS OF BASIS 

FUNCTIONS 

The equation of motion of the covariance were derived by the authors of 0] with the 
final result 

= f "^^C"^ + F"^ , (2.1) 
daC'f = , (2.2) 

where the last equation follows from the solenoidal condition for the magnetic field. 

It is advantageous to to decompose the covariance C"^ in terms of basis functions that 
block-diagonalize the angular part of the operator T. These basis functions are implied by 
the symmetries of T. Since this operator contains only isotropic differential operators and 
contractions with either 5°^ or W^R^, it is invariant to all rotations Q]. Accordingly, the 
natural basis functions should belong to irreducible representation of the SO (3) symmetry 
group, and can be indexed by pairs of indices j, m where j = 0, 1, 2, . . . and —j<m< j. 
We are going to refer to solutions of (|2.1|) that belong to irreducible representation with 
a definite j, m as the "j, m sector". The operator T leaves such sectors invariant. In 
addition, T is invariant to the parity transformation R —R, and to the index permutation 
{a,fi) -v^ {P, u). Accordingly, T can be further block-diagonalized into blocks with definite 
parity and symmetry under permutations. 

In light of these consideration we seek solutions of the form 

C^^iR, t)=J2 a,,UR, t) B-I^{R) (2.3) 

where R = R/R, and -B^j„(/?) are tensor functions on the unit sphere, which belong to 
the sector j,m of the S0(3) symmetry group. The index q enumerates different tensor 
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functions belonging to the same sector. While for scalar functions on the sphere there exist 
only one spherical harmonic Yj^ in each sector, the second rank tensor functions on the 
sphere there exist nine different tensors |Q. The additional symmetries under parity and 
index permutation group them into four subgroups with four, two, two and one tensors 
respectively. With ^jm{R) = R^Yj^iR), in the notation of the 4-group (denoted below 
as subset I) is: 

Bt^R) ^ R-^'-'d'^d^^.UR) ■ (2.4) 

These are all symmetric in a, (3 and have a parity of (—1)-' . The 2-groups are denoted 
respectively as subset II and subset III: 

Bfj^iR) = R-^-\R^e^^^R,d, + R^e''^^R,d,]^j^{R) , (2.5) 
B'^^^^iR) = R-^+^e^^-R^d^d"^ + e''^''R^d,d^]^,m{R) • (2.6) 

B2l^{R) ^ R~^-h"^^R,<^jUR) , (2.7) 
<„(i?) ^ R-^^h-f''^d,^,UR) . (2.8) 

The first couple is symmetric to a, (3 exchange and has parity (—1)-'+^. The second has 
the same parity but is antisymmetric to a,f3 exchange. The remaining basis function is 
B^j^{R) = R~^{R°'d'^ — R^d°')^jm{R) which is antisymmetric to a, (3 exchange, with parity 
(—1)-'. This will be denoted as subset IV. In it was proven that this basis is complete 
and indeed transforms under rotations as required for a j, m sector. 

It should be noted that not all subsets contribute for every value of j. Due to the obvious 
symmetry of the covariance: 

C"^(il,t) =C^"(-i?,t) (2.9) 



it follows that representations symmetric to a, (3 exchange must also have even parity, while 
antisymmetric representations must have odd parity. Accordingly, even j's are associated 
with subsets I and III, and odd j's are associated with subset II. We show below that subset 
IV cannot contribute to this theory due to the solenoidal constraint. 

III. THE MATRIX REPRESENTATION OF THE OPERATOR T 



Having the angular basis functions we seek the representation of the operator T in this 
basis. In such a representation T is a differential operator with respect to \R\ only. In 
Appendix A we demonstrate how T mixes basis functions within a given subset, but not 
between the subsets - as is expected in the last section. In finding the matrix representation 
of T we are aided by the incompressibility constraint. Consider first subset I with 4 basis 
functions Eqs.( |2.4| ) in a given j, m sector. To simplify the notation we will redenote the 
coefficients according to a = agjm,, h = aj^^ c = aijm d = a^jm- In this basis the operator 
T takes on the form 



( \ 

a 




fa") 




1 '\ 

a 




/ \ 

a 


b 




b" 


+ T2 


b' 




b 


c 




c" 




d 




c 














[dj 



(3.1) 



On the RHS we have matrix products. In addition, the solenoidal condition implies the 
following two constrains on a, b, c and d (cf. the Appendix of [§]): 

& , . c 



= a' + 2- + jb' - j^- + c -j- 

X XX 

= 6' + 3- + - + (j - l)d' - (j - - 2)- 

XX X 



(3.2) 



Using these conditions one can bring Ti and T2 to diagonal forms. 
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T2 = -1(M« + k)+?M«] 



T3 can be written in the form 



(3.3) 



(3.4) 



where the four columns of Q{j,i) are 
/ 



(2 + OU + 2)(i + 3) + 2^[{j + 1)(2 + e) + 8] + - 

(2 + 0(2-0 

(2 + 0(2-0(1-0 


^ -2j(j + 1-002-0 

-i(2 + 0(j + i) + 207-0 

-2j02 + 0(2-0 
y 2(2 + 0(2-0 J 

^ -^(2-0(2j-3-0 ^ 

02-0 

-i(2 + 0(i + i) + 0(3 + 
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(3.5) 



-j(j-l)(2-e)(4-0e 
-e(j-l)(2-0(j-4) 

-j(j-i)(2-o(2+oe 

In Appendix B we present the two remaining blocks in the matrix representation of T as a 
function of j. 

Next, the single basis B^jm (subset IV) cannot appear in the theory since asjm = by 
the solenoidal condition (cf. Appendix of [|]): 



3,jm 



jR ^asjrn = 



(3.6) 



Lastly, there are no solutions belonging to the j = 1 sector. This is due to the fact that 
such solutions correspond to subset II. In this subset the j = 1 solenoidal condition implies 
the equation 



''8,1m 



R 







(3.7) 



or ag^im oc R ^ which is not an admissible solution. 



IV. ABSENCE OF DYNAMO EFFECT 

The first issue to clarify is the existence of a stationary solution for t ^ oo. A dynamo 
effect may cause the covariance to grow unboundedly. Vergassola p showed that this is 
not the case in the isotropic sector as long as < 1. We demonstrate that the dynamo 
effect is absent also in the anisotropic sectors. Consider the forceless case of Eg. ( p.lD with 
jpaf3 _ Q addition assume initial conditions such that (B) = 0. It is easy to see that 
no mean magnetic field can appear in time. Accordingly our covariance C"^(i?,t) tends to 
zero when R^ L since C'"^{R,t) {B)"^. We note that for ^ = 0, f"^^ = 2k/\5'^5^. In 
the space of functions C'^^{R,t) which vanish outside the domain \R\ < L this operator 
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is diagonalizable due to its hermiticity, with negative discrete spectrum {—Ex} due to the 
compactness of the domain. Thus the general solution in this case is 



C''^{R,t) = Y.e-^''Cf{R) . 



(4.1) 



In a spherical domain the index A contains the indices j, m and an index specifying one of 
the three subsets discussed above. We will assume that for ^ 7^ T remains diagonalizable. 
We will demonstrate that the eigenvalues E\ remain positive. This will imply that C^^i^R, t) 
is a monotone decreasing function of time, and hence the absence of a dynamo effect. 
To this aim, we define the inner product 



JR 



Ir<l 2{Dm + k) 



(4.2) 



and will demonstrate that 



Ea(Ca,Ca) = fC,,TC>) <0 



A, i ^\) 



(4.3) 



indicating that E\ > 0. We firstly consider the 4x4 block with a given j,m. In this case 
Ca is given by 

Cx{R) = ax{R)Bsj^{R) + bx{R)Bj^^UR) + cx{R)B,j^{R) + dx{R)B^,jm{R) ■ (4.4) 



Using Eg. (^.l]) we obtain 



(Ca,TCa) 



'^dR.rr^ ,(al c*, d^) M{j) 



2{DRi + k) 



K 



( I \ 



yd'xj 



( \ 

ax 



Ca 



ydxj 



where the matrix M{j) arises from the angular integration over the spherical ten- 
sors Bgjrn- This matrix is obtained by a direct calculation. For example Mi i(j) = 
/ dRBlj^{R)Bgjrn{R)- The full matrix reads 
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1 2i 1 j(j-l) 

2j 2j(3j + l) 2j 2j(j-l)(2j + l) 

1 2j 3 



(4.6) 



j{j - 1) 2j(j - l)(2j + 1) j{j - l)(2j - l)(2j + 1) 

We note that M{j) is symmetric and positive definite. By integration by parts, using the 
fact that our covariances vanish for R = L, we demonstrate in Appendix C that Eq. (f4.3|) is 
true. 

One important conclusion of this calculation is the relative rate of decay of the various 
anisotropic contributions. We see that upon increasing j the inner product ( [4.3|) becomes 
more negative. Thus, any anisotropic initial conditions results in a rapid decay of the 
higher j contributions. Without anisotropic forcing the covariance of the magnetic field 
becomes isotropic in time. We will show below that also in the (anisotropic) stationary 
state maintained by anisotropic forcing, the covariance isotropizes on the smaller scales. The 
scaling exponents governing the R dependence are also strictly increasing with increasing j. 
Thus invariably for small enough scales and for long times one restores local isotropy. 



V. CALCULATION OF THE SCALING EXPONENTS 

In the absence of a dynamo effect, we can consider a stationary state of the system, 
maintained by the forcing term f[r,t). The covariance in such a case will bey the following 
equation: 

= f "^^C"^'' + F"^ . (5.1) 

Deep in the inertial range we look for scale invariant solutions, obtained as zero-modes of 
Eq.( [5.1| ). Indeed, when > and well within the inertial range we can take the magnetic 
dissipation to zero, and as a result, the homogeneous part of Eq. ( |5 . 1| ) (without F"'^) will 
be scale invariant, leading to scale invariant solutions. We will need to match these zero 
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modes to the appropriate zero modes computed in the dissipative range at the end. This 
will necessitate the discussion of zero modes when ^ = 0, and see below. 

The calculation of the scale-invariant solutions becomes rather immediate once we know 
the functional form of the operator T in the basis of the angular tensors Bqjrn- Using the 
expansion ( |2.3| ), and the fact that T is block diagonalized by such an expansion, we get a 
set of 2nd order coupled ODE's for each block. To demonstrate this point, consider the four 
dimensional block of T, created by the four basis tensors Sgjm of subset I. According to 
the notation of the last section, we denote the coefficients of these angular tensors in ( p. 3D , 
by the four functions a{R), b{R), c{R), d{R): 



(5.2) 



where . . . stand for terms with other j, m and other symmetries with the same j, m. Ac- 
cording to ( |3.1j ), well within the inertial range, these functions obey: 



Ti(«: = 0) 



a 
h" 



rl" 



+ T2(/« = 0) 



a 
h' 



yd', 



+ Tain = 0) 



a 
b 







(5.3) 



Due to the scale-invariance of these equations, we look for scale-invariant solutions in the 
form: 



i{R) = aR'^, b{R) = bR'^, d{R) = cR'^, d{R) = dR'^ . 



(5.4) 



Substituting ( |5.4| ) into ( [5.3| ) results in a set of four linear homogeneous equations for the 
unknowns a, b,c,d : 

( \ 

a 



[C(C - l)Ti(/€ = 0) + CT2(«: = 0) + T^{k 



. 



(5.5) 
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The last equation admits non-trivial solutions only when 



det [C(C - l)Ti(/€ = 0) + CT2(/€ = 0) + T,{k 



0)] = . 



(5.6) 



This solvability condition allows us to express ( as a. function of j and ^. Using MATHE- 
MATICA we find eight possible values of (, out-of-which, only four are in agreement with 
the solenoidal condition: 



H{^, j) = -e - 8^ + 20'' + 2a + + 4j + 5 . 



Not all of these solutions are physically acceptable because not all of them can be matched 
to the zero mode solutions in the dissipative regime. To see why this is so, consider the 
zero-mode equation deep inside the dissipative regime. Here the dissipation terms become 
dominant and we can neglect all other terms in T. The zero mode equation in this regime 
becomes a simple Laplace equation: 



Notice however, that up to an overall factor of k/{k+D), this equation is identical to the zero- 
mode equation for the special case ^ = 0. Accordingly, the solutions a{R),b{R),c{R),d{R) 
are scale invariant with the exponents Ci''^k=o- Using this we can immediately distinguish 
between the physical and unphysical solutions: At i? = the covariance must be finite, and 
therefore only the positive exponents in ( ^.7|) are acceptable. Let us now consider the zero 
mode equation for ^ = 0. In this particular case, the equation is 



throughout the whole range, and so the zero modes do not change their functional form 
as we move from the dissipative regime to the inertial regime. Therefore for, ^ = 0, also 
in the inertial regime, only the positive exponents are legitimate. Assuming now that the 




(5.7) 



2kV^C = . 



(5.8) 



{2k + 2D)V^C = 



(5.9) 
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solution (including the exponents) are continuous in ^, (and not necessarily analytic!), we 
immediately get, that also for finite ^, only the positive exponents appear in the inertial 
range. Finally there exist two branches of solutions corresponding to the (-) and (+) in the 
square root. 

eg = -\-\^ + yHi^,j) ± 2^Ki^,j) , subset 1 (5.10) 

These exponents are in agreement with Note that in the case of j = only (f^^ exists 
since the other exponent is not admissible, being negative for ^ — 0, and therefore excluded 
by continuity. For j > 2 both solutions are admissible, and the leading one is the negative 
exponent which is smaller. 

In addition to these one needs to compute the exponents corresponding to the subsets II 
and III. The computation in the other two blocks follows the same lines. Since these are 2x2 
they furnish two solutions for the exponents, one of which is negative. We end up finding 

d? = - + + + 2j'e + + 4j + 4j2 , subset II (5.11) 

Q 1 -1 

CS = -- - + 2 V^' + 2^ + 1 + 4j2 + 2f^ + 4j + 2U , subset III (5.12) 

For j = there is no contribution from this subset, as the exponent is negative. The 
dependence of the admissible leading exponents on is displayed in Fig.l and Fig. 2. In 
Table 1 we summarize which are the leading exponents in each sector. 

In addition to matching the zero modes to the dissipative range, one has to guarantee 
matching at the outer scale L. The condition to be fulfilled is that the sum of the zero- 
modes with the inhomogeneous solutions (whose exponents are 2-^) must give C{R) 
as \R\ L. Obviously this means that the forcing must have a projection on any sector 
Bqjm for which a^jm is nonzero. 

VI. SUMMARY AND CONCLUSIONS 

The results of this paper should be examined in light of the recent progress in under- 
standing the effects of anisotropy on the statistics of fully developed turbulence P,|8|-[l0 . 
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Whereas in the Navier-Stokes case one cannot present exact results, the present study can 
affords exact calculations of the whole spectrum of scaling exponents that determines the 
covariance of a vector field in the presence of anisotropy. We have presented a detailed 
and systematic investigation of scaling properties of the covariance of a magnetic field ad- 
vected by a gaussian and delta-correlated in time velocity field. We have extended the 
non-perturbative analysis presented by Vergassola in P] for the isotropic sector to all the 
sectors of the SO (3) symmetry group. Our analysis leads to the conclusions that the scaling 
exponents are strictly increasing with the index of j of the sector, meaning that there is a 
tendency toward isotropization upon decreasing the scales of observation. We also showed 
that as far as the dynamo problem is concerned, anisotropic sectors are less unstable than 
the isotropic sector: in the absence of an external forcing anisotropies decay in time faster 
then isotropic fluctuations. In distinction with the expansion presented in 0, our results 
are free of any assumptions about the hierarchy of scaling exponents belonging to different 
SO (3) sectors. This is due to the employment of a proper basis set. The equations for 
the magnetic covariance foliate into independent closed equations for each set of irreducible 
representations of the SO (3) group. 

In summary, we have shown that the covariance of the magnetic field is naturally com- 
puted as a sum of contributions proportional to the irreducible representations of the SO (3) 
symmetry group. The amplitudes are non-universal, determined by boundary conditions. 
The scaling exponents are universal, forming a discrete, strictly increasing spectrum in- 
dexed by the sectors of the symmetry group. Similar results were presented for passive 
scalar fluctuations in ||^, and for Navier-stokes fluctuations in P,p|-p!0[]. In the present case 
anomalous scaling laws are found as the zero-modes of the inertial operator governing the 
stationary equation for the magnetic covariance PJ^l- Matching with the UV boundary 
conditions selects the physically acceptable solutions. It appears quite clear now that the 
issue of anomalous, universal scaling exponents in turbulence has ramified to the multitude 
of sectors of the appropriate symmetry groups. 



14 



ACKNOWLEDGMENTS 



We are grateful to M. Vergassola for many useful discussions and suggestions. We also 
thank A. Lanotte and A. Mazzino for a fruitful exchange of ideas. At Weizmann this work has 
been supported in part by the German Israeli Foundation, the European Commission under 
the Training and Mobility of Researchers program, and the Naftali and Anna Backenroth- 
Bronicki Fund for Research in Chaos and Complexity. In Rome this research has been 
partially supported by INFM (PRA-TURBO) and by the EU contract FMRX-CT98-0175. 

APPENDIX A: DEMONSTRATION OF THE ACTION OF f"^^ 

As an example of the operation of T on the basis function, consider an explicit calculation 
of d^C"^. Such a term appears as a part of S'^^d'^d" which is a part of T, and also in the 
magnetic dissipation term. Considering explicitly the the part agjrn{R,t)BQ^^{R): 



jm 



d'^a^^jmR-^-^R'^R^^jm = d'^d^agjmR-'-^R^R^^, 

= d>'9JmR~'~' - U + 2)a,JmR-'-']R^^R''R%m 
+d^agjmR~''^R^^jm + 8^ a^jmR'''^ R'^^jm 

- [<,jm - (i + 3)%^ - (j + 2)^ + U + 2){j + 4)^]<^ 
+ (j + 5)[%-(j + 2)^]<^ 

+ ^ U + 2j-^J^9,im + 2-^^l,jm + -^^7jm 

- [<jm + 2% - (J + 2)(J + 3)^]Bf,rn + ^^B^^r. + ^^B^^^ (Al) 

In performing the computation we make use of the following basic identities that are em- 
ployed repeatedly in all our calculations: 

d'^d^^jm = , (A2) 
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(A3) 



The first identity follows from d'^Yjm = + l)R~'^Yjm- The second from the fact that 
$jm are homogeneous polynomials of degree j. As expected, the result remains in a j,m 
sector, and mixes only basis functions with the same symmetry properties. 



APPENDIX B: T AND THE SOLENOIDAL CONDITION IN THE TWO 

REMAINING SUBSETS 

In this Appendix we present the two blocks pertaining to the (—1)-'^^ parity. The part 
denoted in Eq.( p.l| ) as Ti and T2 remains unchanged except that the identity matrix is now 
2-dimensional. For the case of invariance under a, (3 interchange (subset II) we find the 2x2 
matrix Q(j,0: 

\ 



-(j + 1)(2 + i)U + 2 - + 2^(7 - -iU - 1)'(2 - 

(2 - 0(2 + J(J - 1)(2 + + e(j - 3)(2 + + 2^ 

The solenoidal condition reads in this case (cf. Appendix of @]): 

«8jm + 3i?"^a8,jm + (j - l^Jm " (j " ifR'^aQ^j^n = . 



■ (Bl) 



[B2) 



^From this equation we learn that a contribution pertaining to j = 1 cannot appear in this 
theory, since for this value of j a'^j^n must have a negative scaling exponent which is not 
admissible. 

For the case of antisymmetry under a, j3 interchange (subset III) we find the 2x2 matrix 
Q(j,0: 



^ ^(4 + 2^ + 4j) - (j + l)(j + 2)(2 + 0-(j - 1)(2 - ^ 



V 



4-2e -(j-l)[j(2 + + 4e] 

with the solenoidal condition (cf. Appendix of 0]: 



(B3) 



R ^0'A,jrn ~ '^2,jm + (j ~ ^)R ^0,2,jm — 



(B4) 
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APPENDIX C: PROOF OF EQ.(19) 



To demonstrate Eq.(O) we note T as well as M{j) are m independent. We can therefore 
consider without loss of generality the m = case. In this case the basis functions as well 
as the coefficients a, b, c and d are real. For non-zero m the imaginary components have 
to cancel with the imaginary components of — m since the covariance is real. We treat 
separately the contributions associated with Ti, T2 and with T3. showing that they are all 
negative definite. 



1. The integrals of Ti and T2. 

For the evaluation of these integrals it is convenient to work in the basis that diagonalizes 
M{j)- Since M{j) is a real and symmetric matrix, it is diagonalizable and as it is non- 
negative its eigenvalues fii, i = 1, 2, 3, 4, are nonnegative. Ti,T2 are proportional to the 
unit matrix, and therefore they remain so in any basis, and in particular in the diagonal 
basis of M. In that basis, a, b, c, d are replaced by Oi, 02, 03, 04, and the contribution of Ti, T2 
is: 

4^2 / / 

y^,f dx—4 ai [2(Dx^ + K)a'' + A(Dx^ + k)^ + A^Dx^^] (CI) 

~1 J Dx'i + K XX 

This integral is negative definite for all values of i since it is the sum of two negative definite 
integrals Ji and l2'- 
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r d 
-AD / dx-^ 
J dx 





X 



Dx^ + K 



(C3) 



Accordingly, 



-2Dj 



dx 



dx 



X' 



4+1 



Dx^ + K 



= -2D / dx 



Dx^^ + k{1 + i)x^ 
{Dxi + Kf 



at < 



(C4) 



2. The integral of Tg 

The contribution of T3 has two parts: One which is proportional to n, and one which is 
proportional to D. We shall analyze each of them separately and show that M{j) • T3 is a 
nonpositive matrix for every j >2 and every < ^ < 2. 



1. The part involving k is: 



/s = dx 



J ^^ L>J+^ [abed) M{j)KX-^Q{j, 0) 



/ \ 

a 

b 
c 



1j 



d^^^^ {abed) l[M{j)Q{j, 0) + (M(j)Q(j, 0))^ 



/ \ 

a 

b 

c 

ydj 



( \ 

a 



L _^ 



(C5) 



(C6) 



(C7) 



where X (j, ^) is the symmetric matrix 



M(j)Q(j,0 + (M(j)Q(j,0? 



(C8) 



For j = 2, ^ = 0, X{j, ^) is given by 

/ 



X(2,0) = 



-20 -32 -12 
-32 -176 -48 

-12 -48 -36 




(C9) 



with the eigenvalues (—12.97.., —196.43.., —21.59.., 0) so the expression is obviously 
non-positive. For higher j's, we can look at the determinant of X{j, 0): 



det X{j, 0) = (j + 3)(j + 2y{j + irAj - - 2) 



(CIO) 



this function is positive for every j > 2 which means that we have 4 negative eigenvalues 
when j > 2. 



2. The part involving D is: 



/ \ 

a 



L 2 

I-dJ dx^^-^ ( a 6 c d ) M(j>«-2Q(j, 



dx 



x^ 



Dx^ + K 



abed] X{j,C) 



( \ 

a 

b 

c 

ydj 



(Cll) 



The proof of the nonpositivity of this expression follows the same lines of the previous 
discussion. We know that for ^ = 0,j = 2 X{j,^) has three negative eigenvalues 
and one zero. Therefore, It is sufficient to show that detX(j, ^) is positive for every 
< ^ < 2 and j > 2 to ensure that X{j,^) is indeed nonpositive. This is indeed the 
case, as can be verified exphcitly using MATHEMATICA. 
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TABLES 





Symmetric 


Antisymmetric 


J = 


C/+ 




Even j > 


0- 


Cm 


Odd j > 1 







TABLE I. The leading exponents in the various sectors 
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